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Abstract 

We consider the generalized gauge mediated supersymmetry breaking (GMSB) 
models with the messenger fields which do not form the complete multiplets of SU(5) 
GUT symmetry. Such a situation may happen in the anomalous U (1) GUT scenario 
because the mass spectrum of the superheavy particle does not respect SU (5) GUT 
symmetry, although the success of the gauge coupling unification can be explained. 
In this paper, we assume that one pair of the messenger fields gives the dominant 
contribution, and the LHC signature for the two possible messengers, X + X and 
Q + Q, are examined. We investigate the possibility to measure the deviation from 
the usual GUT relation of the gaugino masses which is one of the most important 
features of these scenarios. 
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1 Introduction 



The minimal supersymmetric (SUSY) standard model (MSSM) is one of the most promis- 
ing candidates for the model beyond the standard model (SM) [HEIE]. Unfortunately, the 
MSSM has more than 100 parameters concerning on the SUSY breaking, and the signatures 
in the large hadron collider (LHC) are strongly dependent on the parameters, especially, 
the mass spectrum of SUSY particles [U El CO El El HOI EH]- Therefore, it is important 
to study various possibilities for the SUSY breaking parameters and the LHC signatures 
before the LHC starts to present the data. In this paper, we examine the generalized gauge 
mediated SUSY breaking (GMSB) scenario [121 LT31 EE] in which the messenger fields do 
not form the complete grand unified theory (GUT) multiplets. 

In the usual GMSB scenario [151 [16], the messenger fields are adopted as complete 
multiplets under SU(5). This is mainly because the messenger fields which do not respect 
SU (5) generically spoil the success of the gauge coupling unification in the MSSM. However, 
it has been understood that the GUT with anomalous U{1) gauge symmetry [T71 HHJ fT9] 
can naturally explain the success of the gauge coupling unification in the MSSM although 
the superheavy particles do not respect the SU(5) symmetry [201 EI]- Since some of the 
superheavy particles can play the same role as the messenger fields, it is important to study 
the generalized GMSB scenario with the messenger fields which do not form the complete 
multiplets of SU(5). 

One of the most interesting facts in this scenario is that the GUT relation for the gaug- 
ino masses are spoiled although the anomalous U(l) GUT has the GUT gauge symmetry 
at the GUT scale. Actually, if we assume that the SUSY is broken by some mechanism 
and the resulting superpotential has a form 

W = m*$$ + 9 2 F^, (1.1) 

where 9 2 is a superspace coordinate, the masses of gauginos and sfermions at the scale m$ 
are generated by one and two loop effects of the messenger fields $ and $ as [El [16] 




respectively. Here n a (a = 1,2,3) is the Dynkin index whose normalization is chosen to 
be n a = 1 for 5 + 5 of 577(5), and C[ is the quadratic Casimir invariant of sfermions. C[ 
are in a normalization where C( = 3/5 -Y 2 for sfermions with hypercharge Y , C( = 3/4 

for SU{2)l doublets and C( = 4/3 for SU(3)c triplets. n a for various messenger fields are 
given in Table [1} The generalized messenger scenario has a lot of possibilities in general 
[j~2| [13]. Here, just for simplicity, we assume that one of the messenger fields in Table [T] 
dominates. Then only two possibilities, X + X and Q + Q, can give the non-vanishing 
masses to all the gauginos. In this paper, we examine these two possibilities. 

2 Overview of the mass spectrum 

As discussed in the introduction, the GUT relation for the gaugino masses is generally 
spoiled in the generalized GMSB scenario. Since M a a~ l is one- loop renormalization group 
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Table 1: The Dynkin indices n a (a = 1, 2, 3) for various messenger fields. 



invariant, the gaugino masses satisfy the relation 

Mi : M 2 : M 3 ~ riiai : n 2 a 2 '■ n^a^ (2-1) 

at any renormalization scale. Therefore, if we consider a model with messengers which 
respect SU(5) symmetry, the ratio of gaugino masses at weak scale is given as 

M x {m z ) : M 2 (m z ) : M 3 (m z ) -1:2:6. (2.2) 



The relation (2.2) is often called the GUT relation of the gaugino masses. However, this 
relation is spoiled if ni = n 2 = n 3 is not satisfied as in the generalized GMSB scenario. In 
the followings, we examine the spectra of models with X + X or Q + Q messengers as a 
specific example of such scenarios. The parameters for the SUSY breaking sector are A$ = 
F$/m<s>, m$, tan/3 (the ratio of the VEVs of up-type Higgs and the down-type Higgs), and 
sgn (/x) (the sign of the SUSY Higgs mass). The masses of gauginos and sfermions, and 
the scalar trilinear couplings A at the mass scale of the messenger fields are given as 



M a (m $ ) ~ n a (^) A # , m)(m„) ~ £ n a C[ (^f A|, A(m«) ~ 0. (2.3) 



We use the renormalization group equations (RGEs) to obtain these parameters at the 
weak scale. For implementing the numerical calculation, we use SOFTSUSY 2.0.18 [22] 
with appropriate modification according to our purpose. In our calculation, we assume 
that the contributions from messenger fields other than the selected one are relatively 
small and can be neglected entirely. To specify the SUSY Higgs mass /i and the Higgs 
mixing parameter b we use the relations 



m 2 Hu +m 2 Hd + 2\ii\ 2 )sm2P (2.4) 



and 



l"''= tan^-l -f- (2 ' 5) 

where m 2 H and m 2 H are the SUSY breaking Higgs mass parameters at the weak scale 



which are calculated by the RGEs with the boundary values given by (2.3) at the scale 
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Figure 1: Allowed parameter region for X + X messenger scenario. We set tan /3 = 10 and 
sgn (/i) = +1. The circles correspond to the model points selected to analyze the LHC 
signature. 



First, we consider the scenario with X + X messenger. In this scenario, the Dynkin 
indices are given by 

rii = 5, ri2 = 3, n 3 = 2 (2.6) 
as shown in Table [T], so gaugino masses at the weak scale satisfy the relation 

M x {m z ) : M 2 (m z ) : M 3 (m z ) ~ 5 : 6 : 12 (2.7) 



at one-loop order. Since the relation (2.7) is not affected so much by the specific choice 



of the parameters such as the mass scale of messenger particles, we can use this relation 



to distinguish this model from others. The relation (2.7) indicates that the hierarchy 
among the gaugino masses becomes milder than the usual GUT relation. One of the most 
interesting features is that the mass splitting between the bino and wino is especially small. 
Therefore, to check this feature is one of the promising ways to test this scenario. 

Figure [T] represents the parameter space in this scenario. Here we set tan = 10 and 
sgn (//) = +1, and we assume that the lightest SUSY particle (LSP) is the gravitino. The 
experimental bound for this scenario comes mainly from the LEP2 bound on the lightest 
Higgs mass uih > 114.4 GeV [23]. The mass of Higgs is, however, largely dependent on 
the mass of top quark, so there remains a large uncertainty concerning this bound. We 
set m t = 175 GeV for our calculation. As shown in Figure [TJ there are three parameter 
regions, corresponding to 

1. m^o < m^o < mfj 
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Figure 2: Allowed parameter region for Q + Q messenger scenario. We set tan/3 = 10 
and sgn (/x) = +1. The circle corresponds to the model point selected to analyze the LHC 
signature. 



3. m fl < m^Q < m%o. 

When the A$ is comparatively small, the higgsinos are relatively heavy compared with the 
gauginos, and therefore we can roughly identify the lightest neutralino Xi with the bino 
and the second lightest neutralino x® with the wino. Therefore, the bino-like neutralino 
Xi becomes the next to LSP (NLSP) in the regions 1 and 2, and the stau T\ becomes the 
NLSP in the region 3. 

Another candidate of messenger fields in our scenario is the fields with quantum number 
of Q + Q. In this scenario, the Dynkin indices are given as 

m = -, n 2 = 3, n 3 = 2 (2.8) 
5 

and gaugino masses satisfy the relation 

M^mz) : M 2 {m z ) : M 3 (m z ) ~ 1 : 6 : 12 (2.9) 

o 

at one- loop order. As can be seen from this relation, Q + Q messenger scenario gives rather 
small masses to sparticles which do not have quantum numbers of SU{2) L and SU(3)c 
compared with the other sparticles. The allowed parameter space is presented in Figure 
|2j In the whole allowed region, m^o < < m^o is satisfied as in the region 2 of X + X 
messenger scenario, and therefore the NLSP is the bino-like neutralino Xi- 
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Note that there are almost no bounds for the mass of the lightest neutralino Xi if Xi 
is the pure bino and does not decay inside the detector. The GUT relation is essential to 
obtain the bound m^o > 46 GeV given by the particle data group [23] and the constraint 
from the invisible decay of Z is useless because the decay width of Z — > Xi Xi is quite 
small for the bino-like Xi [21] • Therefore the constraint for the mass of right-handed stau 
fi, m f > 81.9 GeV is important in this scenario which is shown in Figure [2j 



3 LHC signature 

In this section, we investigate the LHC signatures of these scenarios. For this purpose, 
we use ISAJET 7.79 [25] to calculate the decay width of sparticles and HERWIG 6.510 
[26] 127] to generate the sparticle production events by Monte-Carlo simulation. And for 
the detector simulation, we use AcerDET 1.0 [2H] as a fast simulation of the search at the 
LHC. We examine the LHC signatures for y/s = 14TeV for the whole analysis in this 
paper. 

We pick three model points for X + X messenger scenario corresponding to the three 
regions introduced above (Table [2J and one model point for Q + Q messenger scenario 
(Table [3]) to analyze the LHC signals. Table [2] and [3] show the resulting mass spectra and 
branching ratios of sparticles in these model points. In these points, the NLSP does not 
decay to the LSP gravitino inside the detector. 

As pointed out above, one of the most peculiar features of these scenarios can be tested 
by measuring the masses of the bino and wino. For the X + X messenger scenario, the 
mass splitting between the bino and wino is very small. On the other hand, for the Q + 
Q messenger scenario, the mass of bino is much smaller than other sparticle masses. So 
one of the most important tasks to distinguish these scenarios is measuring the neutralino 
masses 

m^o ~ Mi, m^o ~ M 2 . (3.1) 

Of course we have to measure the mass of gluino to confirm the relation among gaugino 
masses predicted by our scenarios. But we do not argue the detailed reconstruction of the 
decay chain for the gluino mass measurement in this paper. This is because the decay 
mode of the gluino in the GMSB model is highly dependent on the parameters and it 
can be very complicated. As shown later, the gluino mass can be estimated if we assume 
that the gluino mass is of the same order as the squark masses, for example, by the rriT2 
measurement and the largeness of the cross section. 

We figure out several features for three cases in X + X messenger scenario and one case 
in Q + Q messenger scenario. 

1. : irig ~ 1 : 2 in both scenarios. As the result, the ratio m^/rriq becomes larger 
than in the usual scenario with the GUT relation. Roughly speaking, the hierarchy 
between colored sparticles and wino masses becomes milder. 

2. nig ~ m^r in the X + X messenger scenario. In most of the interesting parameter 
region, the decay mode x\ ~~ >" %Xi is closed, and therefore, the branching ratios of 
leptonic decay modes of X2 become comparatively large. The leptonic modes are 
important in obtaining meaningful information from the data in the LHC. 
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Table 2: Mass spectra and branching ratios of sparticles for three model points correspond- 
ing to three regions of X + X messenger scenario. (See Figure [TJ) 



6 



-A-mess 


^mess 


tan (3 sgn (/i) 


n\ n 2 n 3 


(TeV) 


(GeV) 






80 


1.0 x 10 14 


10 + 


1/5 3 2 



g 1181 Xi 16 u L ^~gu .25 

u L 1326 xl 586 ^ xt d .39 

Wtf 1165 Xs 682 ^xlu -19 

J L 1331 xl 720 ~d~^j~d ^6 

rffl H63 Xi 586 u 35 

§i 1150 *2 720 ->x°rf .18 

6 2 1221 /i° 116 7^->Xi« 1-00 



f 2 1240 A 954 
z>, 680 i? ± 958 
i> T 679 



1.00 



g^ht + tlt .97 

r,o v ,~o 



R8/I X2->Xi^ -30 

e ^ 684 Ai ^1 

e~ R 118 ~+4 h - 54 

f, 96 ->^ /T 

f 2 682 l^X?^ 1-00 



Table 3: Mass spectrum and branching ratios of sparticles for a model point of Q + Q 
messenger scenario. (See Figure [2j) 



3. <C in the Q + Q messengner scenario. 
In the followings, we study how to catch these features from the LHC signals. 



3.1 X + X messenger scenario (Case 3: stau NLSP) 

In the region 3 of X + X messenger scenario, very peculiar signal is expected because the 
NLSP becomes the right-handed stau. The momentum and velocity of stau which goes 
out the detector can be measured, and therefore, we can know the masses of various SUSY 
particles by using the invariant mass technique and the above features can be tested. We 
discuss the mass measurements for this case in this subsection. 

It has been studies how to catch the stau in the LHC in (32j |33j Ell EH], and we carry 
out the smearing of the stau momentum and velocity to reproduce the expected resolution 
in our simulation. The resolutions for the momentum and velocity are given as 

o-|p|(GeV) = 0.000118 • \p\ 2 + 0.0002 • ^jm\ + |p| 2 + 0.89 (3.2) 

and 

a/? = 0.028 • f3 2 . (3.3) 

Then we can obtain the stau mass 

m fl = |i, 7 = (l-/3 2 )-V2 (3.4) 
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Figure 3: Left: NLSP stau invariant mass distribution for the case 3 of X + X messenger 
scenario. Center: m?i distribution. Right: Hemisphere invariant mass distribution. 



from the reconstructed momentum and velocity. For the identification of stau, we use the 
cuts 



0.9 < M < 6-0 
\r]\ < 2.4 
p T > 20 GeV 



in our analysis. Here 77 is the pseudo rapidity. 

To measure the neutralino masses m^o , we use the decay chain of Xi 2 — ^R ~~ >" 
f\ l n lf t. Since mass difference between the lighter stau f\ and the first two generation of 
right-handed sleptons Ir are small in most of the cases we are interested in, we use the 
following approximation: 

PlR =Pr+ Plf +Pr~ Pfl ■ ( 3 - 5 ) 

So we can measure the masses of the bino-like and wino-like neutralinos directly by taking 
the invariant mass of stau and e or /i as illustrated in Figure [3] Although we can use the 
decay chain X12 ~^ f° r the mass measurement of neutralinos, by this measurement the 
events with e or fi are preferable because the momentum of e or /1 is less smeared than 
that of r. And we can confirm that the measured neutralinos are not higgsino-like one 
because of the large branching ratio for X12 ~* ^R^- 

In this scenario, it is expected that we can measure the masses of any kind of sparticles 
produced in each event. This is because the momentum of the arbitrary sparticle can be 
reconstructed by the momenta of the stau and the SM particles. So we can measure the 
mass of gluino and check the mass relation among all the gauginos in principle. Let us 
consider this issue in the rest of this subsection. 

Although a large number of gluino are expected to be produced by the process pp — > 
gg, gq and the subsequent decay of squark, it is not always easy to distinguish gluino 
from squark in event-by-event level. Therefore we adopt the inclusive measurement of the 
invariant mass of produced sparticles. For this purpose, we use the hemisphere method 
suggested in [HI EHJ [37] . In this method, we sort the clusters into two hemispheres in each 
event according to the following algorithm. 
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1. We pick all the jets with px > 50 GeV, \rj\ < 2, leptons with p T > 15 GeV, \r]\ < 
3 and two staus in each event. We use them as the clusters which compose two 
hemispheres corresponding to the pair-produced sparticles. 

2. We define the initial hemisphere axes p^ erai (i = 1,2) by the momentum of two 
clusters. p£2ni * s defined as the momentum of the highest pr cluster. And PheLi 
corresponds to the momentum of the cluster which has the largest value of pt x AR 
where AR = a/ ' Arf + A0 2 . Here Arj = rf^ — rf 2 \ Acj) = <fP> — <j)^ and is the 
azimuthal angle of the cluster. 

3. The cluster with momentum p is belonging to the hemisphere 1 if it satisfies 

<*(P>Ph2ni) < ^PheLi) ( 3 - 6 ) 

and vise versa. Here d(pj,p^l mi ) is the Lund distance measure between the clusters 
with momentum pj and p^l mi , and it is defined by 

d(Pi,p£L) = (4L - lp£Ll cos%) -J^—- (3-7) 

where 0^ is the angle between and Phemi- 

4. We redefine the hemisphere axis p^mi as ^ ne sum °^ ^ ne momenta of the clusters 
which belong to the hemisphere i. 

5. We repeat the step 3 and 4 until the classification of hemisphere converges. 

After this algorithm, we can obtain the invariant mass distribution of each hemisphere 
m hemi = Phemi- ^ ^ e assignment of hemisphere agrees with the true hemisphere, mh em i 
corresponds to the mass of the pair produced sparticle. Since the true hemisphere should 
contain exactly one stau, we reject the event where two staus are contained in one of the 
hemispheres. 

Using this algorithm, we illustrate the distribution of mh em i in Figure |3j We can see 
from this figure that both gluino and squark have masses around 1000 GeV and the mass 
relation : m~ g ~ 1 : 2 can be checked. 



3.2 Low luminosity analysis 

In order to catch the features in the other cases of X + X messenger scenario and of Q + 
Q messenger scenario, we discuss the several analyses which can be done in the LHC with 
comparatively low luminosity. 

Since we are interested in the measurement related to the masses of the bino and wino, 
we make use of the characteristic decay modes of these particles. As noted in the comment 
of feature 2, for X + X messenger scenario, the leptonic decay of wino-like neutralino x° to 
the lightest bino-like neutralino Xi through on-shell or off-shell slepton Ir is useful, which 
are illustrated in Figure |4| 

First, we consider the invariant mass of a pair of leptons coming from the decay of 
X2- As can be seen from Table [21 in the X + X messenger scenario, a large number of 
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Case 1: / I Case 2: 




Figure 4: Decay chains which we use for the analysis in this paper. X2 decays to Xi y i a 
the off-shell slepton in the case 1. In the case 2, X2 decays to on-shell slepton and we label 
the leptons which come from the decay of X2 an d Ir as l n and //, respectively 
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Figure 5: Dilepton invariant mass mu distributions for 10fb _1 of each model points. The 
bins illustrated by the dotted line represents the similar distribution using the opposite-sign 
and different flavor dileptons (e + fi~ + /i + e~). 



X2 is expected to be produced via the decay of left-handed squark q^. We use a dilepton 
with opposite-sign and same-flavor (e + e~ + fi + fi~) for making the invariant mass and the 
distribution shown in Figure [5] is obtained. In this figure, the background is estimated 
by taking the similar invariant mass of dilepton with opposite-sign and different-flavor 
(e + /i~ + /i + e~). Here, in order to reduce the SM background, we impose following event 
cuts by using the transverse momenta pr [5] 

• p { f > 100 GeV and ^ 2,3 ' 4) > 50 GeV 

• M eff = + p { p + pf + pP + Ef ss > 400 GeV 

• E™ iss > max{100GeV,0.2M eff } 

• Two isolated leptons with p e T > 20 GeV and p^. > 5 GeV 

where pj) means the i-th largest px of the jet in each event and E™ 1SS = -^/(p™ lss ) 2 + (p™ iss ) 2 - 
Since the SM background is reduced successfully after these cut, we generate only events 
of sparticle production for our simulation [5]. 

We can see the rather small maximum value of invariant mass for both cases of the X + 
X scenario in Figure [5J which is caused by the feature 2, namely, ~ m-w- Actually, 
the maximum value of invariant mass allowed by kinematics is given as 

m™ ax [Case 1] = m^o - m^o (3.8) 

in region 1, and 
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Figure 6: px distributions of leptons emitted from sparticle decays. Each graph corresponds 
to lfb" 1 for the case 1 of X + X scenario (left), lfb -1 for the case 2 of X + X scenario 
(center) and 10 fb -1 for Q + Q scenario (right). 




'»<;r[C*,se2], 2 I 1 - T3 1 ) '3.9) 

in region 2, which result in rather small maximum value of the invariant mass calculated 
as 

m™ ax [Case 1] = 44 GeV, m™ ax [Case 2] = 47 GeV. (3.10) 

Unfortunately, the smallness of m™ x does not always mean the smallness of the mass 
splitting between the bino and wino. For the case 2, if one of the relations ~ m z ~ and 
mj ~ rriyy is satisfied, the maximum value of the invariant mass becomes small. 

To check that the mass difference between \\ an d X® is small in X + X scenario, we 
examine the px distribution of leptons which come from the decay shown in Figure |4| Since 
the magnitude of pt of produced particles is strongly dependent on the mass difference 
among the sparticles, this can be a good signal to distinguish these scenarios. We can see 
from Figure [6] that these leptons have relatively small pt in the X + X scenario and large 
Pt in the Q + Q scenario. These are caused by ~ in the X + X scenario and 
mg <C rriyy in the Q + Q scenario. 

Even in the models with the GUT relation for the gaugino masses, such a small value 
of m™ ax is possible if the gaugino mass scale is small. But such lighter gluino can be distin- 
guished from heavier gluino in X + X messenger scenario by measuring the cross section 
and/or by the m^2 method. As a reference model, we adopt the minimal supergravity 
(mSUGRA) model with parameters mx/ 2 = 150 GeV, m = 750 GeV, A = —100 GeV, 
tan/3 = 10 and sgn (/i) = +1. As shown later, though the distributions of mu and px in 
the reference model are similar to those of X + X messenger scenario in Figure [5] and |6j 
the distribution of rriT2 and the cross section become much different from those of X + X 
messenger scenario. Let us remind the 7717-2 method [291 EH EI] . When we consider the 
production process of sparticle pair which decay into a pair of the NLSPs and a pair of 
the SM particles, we can make use of the 777,7-2 variable defined by 

777 T2 (M test ) = min max \ m T Cf^ sil \ ~?^ (1) ), m T (~^! s{2) , ~f^ (2) ) \ . (3.11) 

pmiss =Ep x(0 L I J J 

i 

Here ~$t is a (2 + l)-dimensional vector, ~$t = (Et,Pt), and p™^ is the transverse 
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Figure 7: ra T2 distributions for M test = 0. 
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momentum for the emitted observed particles, my is given by 



ran 



<J^\ -fe) ee (-ff + ^ T f = m*. + M t 2 est + 2(£™ • - p* ■ p*) (3.12) 



where E™ 



rat 



|p™| 2 and E% = ^/M^+JpW- 



.(2) 



Now we consider the process pp — > q q — > Xi q xi Q- In this process, ~ ra^ n ~ 



and the maximum value of ra^ is given by 



ra\ - m| 
m£ 2 ax (M tcst ) = 9 2m _ Xl + 



\ 



2 2 \ 2 

ra- — m~ Q 
" xl 1 + M? 



2m,- 



tcst 



(3.13) 



as a function of the test mass M tes t- Therefore, we can obtain the rough value of the 
colored sparticle masses by this analysis. 

For the analysis of rriT2i we use ^ ne event cuts 

• Two jets with p T > 100 GeV 

• M eff = p^ + pf + Ef ss > 400 GeV 

• E™ iss > max{100GeV,0.2M eff } 

• No lepton 



instead of the usual cut for the SM events introduced above. Since (3.13) is satisfied for 
any fixed value of M tes t, here we set M test = for our analysis. Then (3.13) becomes the 
following simple form: 



ra- 



m1 

Xi 



rag 



(3.14) 



for the process of squark pair production. The distribution of this quantity is shown in 
Figure [7| If the mass of the LSP is very small compared with the mass of squark, we can 
interpret m™2 X (0) as the mass scale of squark. In fact, this is the case for Q + Q messenger 
scenario. However, the mass hierarchy of sparticles is small in X + X messenger scenario 
and the effect of ra^o is non- negligible. The theoretical values of m™| x (0) are 



m 



T2 



: (0)[Case 1] = 844 GeV, 



m 



T2 



: (0)[Case 2] = 904 GeV, 



(3.15) 
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Figure 8: The measurements for the mSUGRA model point, where m^j = 150 GeV, m = 
750 GeV, A = —100 GeV, tan (3 = 10 and sgn (/i) = +1. Left: Dilepton invariant mass mu 
distribution. The bins illustrated by the dotted line represents the similar distribution using 
the opposite-sign and different flavor dileptons (e + + fi + e - ). Center: px distribution of 
leptons emitted from sparticle decays. Right: m^2 distribution for M tes t = 0. 



here we approximate rrig = 1000 GeV for the case 1 and rriq = 1100 GeV for the case 2. 
By this analysis, we can obtain the evidence of milder hierarchy between the masses of X2 
and colored sparticle if we know m^o ~ 500 GeV. Unfortunately, we do not find the scale 
of m^o by the analysis in this subsection and it needs further detailed analysis. 

In the reference model, the distributions mu, pt and m^2 are represented in Figure |8j 
The distribution of m^i is much different from those in Figure [7[ although the distribution 
of mu and pt is similar to those of X + X messenger scenario in Figure [5] and |6j And the 
cross section becomes much larger than in X + X scenario. Since the distribution of 
and the cross section in X + X messenger scenario show the much larger mass scale of the 
colored particle than the gluino mass obtained by the GUT relation, it is suggested that 
the GUT relation is not satisfied. 



3.3 X + X messenger scenario (Case 1, 2: neutralino NLSP) 

In the cases 1 and 2 of X + X messenger scenario, the neutralino is the NLSP and we 
focus on the decay chain shown in Figure |4j Although these two cases give similar signals, 
the decay modes of the wino-like neutralino are different. In the case 1, x® undergoes 
three-body decay through off-shell slepton, while x° decays to the right-handed slepton 
Iji, which decays to Xi subsequently, in the case 2. As seen in the previous section, the 
small values of m™ ax and p l T indicate that m^o — m^o is small, but this may not mean 
that m^o ~ m^o because the possibility may be still alive that the absolute value of the 
neutralino mass scale is small. In order to reject the possibility, we try to show the relation 
m^o ~ m^o by measuring the invariant mass mji^ of a jet emitted from the squark q and 
one of two leptons in the decay of X21 although the large luminosity is required for this 
analysis. Since there are two leptons in each event, we include two invariant masses mji^ 
for each event in the distribution. In the case 1, the maximum value of mji( u \ is obtained 
as 

MS [Case l]) 2 = ml ( 1 - ^§ ) ( 1 - ^ ) , (3.1 6 ) 
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which is predicted to be 392 GeV for rriq = 1000 GeV. Note that this predicted value is 
much smaller than m™2 X (0) ~ 844 GeV. This indicates that m^o ~ m^o unless rriq ~ r^o- 
Note that in the case 2, there are two kinds of leptons in the decay because there is an 
on-shell slepton Ir produced by the decay of Here we label two leptons emitted from 
X2 and Ir as "near"-lepton /„ and "far"-lepton If, respectively, as shown in Fig 
we cannot distinguish l n with If in event-by-event level, we consider the quantity 



4j Since 

mji^ = m fln U m flf (3.17) 



suggested by [3"""", 159*] . mji n {mji f ) means the invariant mass of a jet emitted from squark 
q and a lepton l n (lf). Then rriji^) gives a combined distribution of rriji n and rriji r The 
important point is that the analysis is completely the same as in the case 1, and we do 
not have to distinguish l n and If at event-by-event level. The maximum value of the 
distribution of Triju u ) becomes 

max{mX x ,m- ax }, (3.18) 

where 

KT) 2 = < I 1 - dr I I 1 " dr I (3-19) 

and 

KT) 2 = m '~ I 1 - df ) I 1 - dr I • (3-20) 

In the case 2, is predicted to be 369 GeV because m™ ax = 262 GeV and m™ ax = 369 

GeV for rriq = 1150 GeV. Again, this predicted value is much smaller than m™2 X (0) ~ 904 
GeV, and it indicates that m^o ~ m^o because rri 1 ^ <C m™2 X (0) means m^a ~ mj- and 
m ™^ x ^ ^t| x (0) m eans m^o ~ unless rriq ~ "T.^' Therefore, if Tnf^, is much smaller 
than m™| x (0), m^,o ~ m^o can be shown. 

In our simulation, we impose the cuts for the standard model background as in the 




section 3^ and use the dilepton whose invariant mass mu is less than 50 GeV and a jet 
with pt larger than 100 GeV. There are, however, many background of jets coming from 
other decays of colored sparticles, such as g — > t\ t — > & t- To reduce these background, 
we impose another event cut that there is no 6-tagged jet in each event. Here we assume 
60% tagging efficiency of 6-jet. Then we make two invariant masses (i = 1,2) for 
all the possible jets in each event. And we choose a jet which minimizes max{rriji 1 , rriji 2 } 
among these jets. In this way we can obtain the distribution of w,ju u \ which consists of the 
combined distribution of and rriji 2 . The above predicted values are roughly consistent 
with the measured values obtained from Figure 19} 



Note that in the case 1, taking account of three relations (3.8), (3.14) and (3.16) 
together, we can obtain the masses of squark q and neutralinos x° 2 , in principle. We will 
return to this point in the end of this subsection. 

In order to determine the mass spectrum of sparticles in the case 2, more measurements 
are needed. As suggested by [39], min{m^ x ,m^ x } can be obtained from the distribution 
of mju u ) by observing the structure with the intermediate endpoints in the case 2. Actually 
the predicted value roughly agrees with the measured value in Figure [9] 

Although we do not know which of the two measured endpoints of rrijiu^ in the case 
2 corresponds to rriji n (mji f ), it is shown in |39j that we can determine the masses of 
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Figure 9: m,ji( u ) distributions for 300 fb of the case 1 (left) and case 2 (right) of X + X 
scenario. 



neutralinos without ambiguity by introducing the quantity 



The maximum value of this quantity is obtained as 



(m^[Case 2]f 



2 

m~ 








+ R 




) 2 










m\ J 







= m|| l--f )( 1-^). (3.22) 

in the case 2. Then we impose a similar event selection as in the previous one and make 
an invariant mass by a dilepton with mu < 50 GeV and a jet with px > 100 GeV. Among 
all the possible choice of a jet, we pick the one which minimizes m^( s ) in each events. And 
the result of simulation is shown in Figure 10, whereas the theoretical value is calculated 

as 

mf^ [Case 2] = 442 GeV. (3.23) 

From these quantities, we can obtain the neutralino masses without any ambiguity in the 
case 2 [39]. By denoting that 

A = m™ ax , B = max{mJ^,m^}, C = min{m^ x , mf^}, D = m^, (3.24) 

the masses of neutralinos are written as 

i Aa/ (D 2 — B 2 )(D 2 — C 2 ) r ^ . ABC 

m^o[Case 2] = B 2 + C 2 _ D 2 > m^o[Case 2] = B2 + C2 _ D2 - (3.25) 

Moreover, the squark mass is also obtained as 

ffCyM 2 + £? 2 + C 2 — £> 2 
m 9 -[Case 2] = ^ 2 - g2 _ ^ . (3.26) 

If we take A = 50 GeV, 5 = 380 GeV, C = 240 GeV and D = 440 GeV, we obtain rn^o ~ 
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Figure 10: rriji^ distributions for 300 ftT 1 of the case 1 (left) and case 2 (right) of X + X 
scenario. 



490 GeV, m^o ~ 540 GeV and m% ~ 1130 GeV, which are in good agreement with the real 
values in Table [3] The measured value of m™| x (0) can be used to check the consistency. 

We can also calculate the maximum value of nijUs) for the case 1 as 



m 2 ~ 
2 i xi 



m^[Case 1] = m™ x [Case 1] = mi 

= m^ ) [Ca 8 el]. (3.27) 




rfus) = m ll - m l ( 3 - 28 ) 



where rriju means the invariant mass by the jet and dilepton. The first equality can be 
shown by the trivial relation 

"'Ms) = " l W 

and the fact that rriju is maximized when m\i = as long as m^o/m^o < rriq/m^o is satisfied 
[38] ■ Let us determine the masses of squark and neutralinos Xi2- For example, if we take 
m max = 45 GeV, m^ 2 ax (0) = 800 GeV and mf^ = = 400 GeV, then we can obtain 

m^o ~ 330 GeV, m^o ~ 370 GeV and rriq ~ 920 GeV which are not far away from the real 
values in Tabled 



3.4 Q + Q messenger scenario 

In Q + Q messenger scenario, we can see the dilepton signal if we collect a large number 
of events. The maximum value of m« is given as 

1 - (l - „f J (3^29) 

and the predicted value becomes 

m^ ax = 569GeV (3.30) 
in this model point. The result of simulation corresponding to 100fb _1 is shown in Figure 



11 and the measured value is consistent with the predicted value. Since m™ ax gives the 
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Figure 11: Left: Dilepton invariant mass mu distribution for 100fb _1 of Q + Q messenger 
scenario. The bins illustrated by the dotted line represents the similar distribution using 
the opposite-sign and different flavor dileptons (e + /i~ + ji + e - ). Right: leptonic m T2 
distribution for 100fb _1 of Q + Q messenger scenario. 



lower bound of the mass of X21 such a large value of m™ ax indicates the large m^o. In 
the section |3~2 we have already the mass scale of the heaviest colored particle, which is 



predicted as 1200 GeV, by analysis for the colored sparticle pair production. These 
signals mean that the mass ratio m^o/rrig is roughly larger than 1/2, and therefore, the 
hierarchy between the gluino and wino masses is milder than that in the models with GUT 
relation. 

We can also use the mr2 analysis for the right-handed slepton pair production, because 
the right-handed slepton also has a smaller mass compared with the other sparticles. Using 
the two leptons emitted from a pair of sleptons, we obtain the maximum value of m<r 2 
variable for M test = GeV as 

,,, m~ — m 2 - 
m « max (0) = *i (3.31) 

Ir 

To select the process pp — > — > Xi l + Xi l~ > we impose event cuts so that 

• No jet with p T > 20 GeV 

• Two leptons with p 1 ^ > 50 GeV and no other leptons with p l T > 20 GeV 

• Ef ss > 50 GeV 

• p f^ + p 1 ^ + £;miss > 200 GeV 

• The invariant mass of two leptons is outside the region 80 GeV < mu < 100 GeV 

are satisfied. Note that the last cut is imposed to suppress the SM background. Since 
the expected background coming from the SM events is leptonic decay of XV boson and 
Z boson, we have to consider the m^ 2 distribution for WW, ZZ, and ZW production 
processes. For the process pp WW —> lulu, it is obvious that 
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mj 2 (0) < m w (3.32) 

and the process pp -^ZZ^llvvis rejected by the above cut. For pp — >■ ZW — )■ lllv ', 
this process will contribute to the m-T2 distribution if one of the leptons comes from Z 
boson is missed to be detected. In that case, 

m^O) < max{mz,mw} = mz- (3.33) 

Therefore, the standard model background does not affect the measurement of m^ m!K (0) 
if (m~ — m1 )/mr > m? is satisfied. We checked that these SM background can be 
negligible by producing the SM background from the ZZ, WW, and ZW processes with 
the above cuts. 

In the model point we are considering here, the maximum value of this rriT2 variable is 
given as 

M = m? 2 max (0) = 116 GeV (3.34) 



and the result is shown in Figure 11 Here we illustrate the distribution in the region 
m^ 2 (0) > m z , because the SM background is expected to be negligible only in this region. 

The important relation m^o 3> m^o can be obtained by the calculation as 

m-o m 2 ~ m 2 ~ — m~ (m maxS \ 2 

_j£ > _*t > _*1 in = l™« ) „ 2 g TeV; (3 35) 

m x? m i R m i R m^ 2 max (0) 

where we use the measured values of mjj 121 * and m^ 2 max (0). 

It may be possible to measure the masses of mj R and m x o by the methods discussed in 
the papers [HI H2] • Here we do not discuss this issue further. 



4 Discussion 

Generically, multiple fields may play as the messenger fields. Actually, any vector-like 
fields $j and $j (i = 1, • • • , n) can be the messenger fields if they have an interaction with 
the spurion field S like n^^iS. Here the F component of S has non-vanishing VEV F$ 
which breaks the SUSY. Naively, they have the same order of the contribution to the SUSY 
breaking parameters if the coefficients Ki ~ Qm^JA where A is the cutoff. Namely, the 
scale A$ = F$ i /m$ i ~ a^Fs/A becomes independent of i except the coefficients q of the 
interactions m$ i $j$jS'/A. Because of the freedom of the 0(1) coefficients, we have various 
possibilities for the sparticle spectrum. Indeed, any gaugino mass spectrum are possible 
by choosing the coefficients Cj. In this paper, we consider an extreme case, in which one 
of the coefficients becomes much larger than the others. This can often happen when the 
coefficients are the ratio of the 0(1) coefficients q = at/bi, where the 0(1) coefficients eij 
and bi are determined randomly, for example, between and 1. Since it is reasonable to 
expect that one of the n coefficients of the denominator, which is noted as b\, becomes 
0(l/n), the coefficient C\ can be 0{n). If n ^> 1, very large coefficient c\ is realized. 
Moreover, if b\ happens to be 0(l/(lOn)), which requires 10% tuning, c\ ~ 0(lOn), and 

n 

therefore, the messenger field $1 and $1 can dominate the others. If °i ~ 0(100), then 

i 

since the loop suppression is almost compensated by the summation of the coefficients, the 
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contributions to the SUSY breaking parameters from the gauge mediation may become of 
the same order as those from the direct interactions between the spurion field and MSSM 
fields, for example, Q^QS^S/A 2 in the Kahler potential, n can be much larger than 1, for 
example, in the anomalous U(l) GUT scenario, there are many vector-like fields which 
have non-trivial charges under the standard gauge group and can be the messenger fields. 
There are 23 vector-like pairs in the S'O(IO) model [T?] |2"U], 47 pairs in the E Q model [T£] , 
38 pairs in the simpler Eq model [T9] except the MSSM Higgs doublet pair. 

In this paper, we consider the several cases in which one vector-like messenger field 
dominates the others. Note that although we have chosen X+X or Q + Q as the messenger 
field in order to obtain non- vanishing gaugino masses, other choices become possible if the 
small but non-zero contributions from the other vector-like fields are taken into account. 

If there are multiple messengers $j and (i = 1, • • • , n) which has the same quan- 
tum numbers, D-term contribution of U(l)y hypercharge to the sfermion masses becomes 
important [10]. This contribution comes from one-loop Feynman graph and may not be 
negligible. The explicit formula is given by 



12 

f 2 z ^\AnJ f * j-f max{(m$J 2 , (m^.) 2 } 

where F$y is a SUSY breaking mass mixing parameters (so called B parameter) of the 
messenger fields in the unit in which the mass matrix of the messenger fields is diagonalized 

as 

^ mess = m^A + <IyI' r (4.2) 

If the enhancement factor c 3> O(10 2 ), then such contribution can be negligible. 



5 Summary 

In this paper, we investigated the LHC signatures of the generalized GMSB models with 
the messenger fields which do not respect SU(5) GUT symmetry. Such a situation can be 
realized in the anomalous U(l) GUTs in which the success of the gauge coupling unifica- 
tion can be explained although the mass spectrum of the vector-like fields do not respect 
SU (5) GUT symmetry. The mass spectrum of sparticles become different from those in 
the usual GMSB whose messenger fields respect SU (5) GUT symmetry. Especially, the 
gaugino masses do not satisfy the usual GUT relation and this feature is very important 
to distinguish these models by the LHC measurements. In principle, any mass pattern for 
the gaugino masses is possible in this generalized GMSB scenario. In this paper, only for 
simplicity, we examined the models with a pair of messenger fields which have quantum 
numbers of X + X or Q + Q and studied how to obtain the signatures of these models in 
the LHC. The gaugino mass relation becomes : : rrig ~ 5 : 6 : 12 for the X + X 
messenger model and : : rrtg ~ 1/5 : 6 : 12 for the Q + Q messenger model. One 
of the interesting features of the both models is that the hierarchy between the colored 
particle masses and weakly charged particle masses becomes milder than the usual GMSB 
models because the messenger fields have bi-fundamental representation under SU(3)c x 
SU(2)l. If we catch the scale of the wino mass in the LHC, we can roughly check this 
milder hierarchy by measuring by which the order of the colored particle masses can 
be obtained. In X + X messenger scenario, the mass hierarchy between the bino and 
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wino is very small, and it leads to the relatively soft px distribution of leptons. On the 
other hand, the mass difference between the bino and wino is considerably large in Q + 
Q messenger scenario. Therefore very large pt distribution of leptons can be seen at the 
LHC. 

The NLSP of X + X messenger scenario is the stau or neutralino. If the NLSP is 
the stau, we can check the mass relation of gauginos at the low-luminosity stage of the 
LHC and the deviation from the usual GUT relation can be obvious. If the NLSP is the 
neutralino, we can determine the bino-like and wino-like neutralino masses by the use of 
the end-point in the neutralino's leptonic decay and of the measurement. 

In Q + Q messenger scenario, the leptonic m^2 measurement is useful because the 
right-handed slepton remains light. 

Since both scenarios predict very characteristic mass spectra, it is expected that we 
can distinguish these models from the models which satisfy the GUT relation. 
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A Two-loop RGE effects for the gaugino mass 



There is a well-known feature in the softly broken SUSY models. Namely, by looking 
one-loop RGEs of gauge couplings and gaugino masses, one can see that their ratio obeys 
the following RGE. 

Then, the gaugino masses in the GMSB model satisfy the relation 



: AI-, : M 3 



ni«i : n 2 a 2 : n 3 a 3 



(A.2) 



at any renormalization scale. But it should be noticed that (A.l) is satisfied only up to 



one-loop order, so the actual relation is deviated from (A.2) by two-loop order correction. 



Let us estimate the effect of two-loop order RGE. The two-loop RGE of gauge couplings 
and gaugino masses are given by 



d 



ft, 



d\nQ a 
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'2vr 



EB a b c a 

b Stt z 



(A.3) 



and 



d 



dlnQ 



-Ma = ^-a a M a + 

2n 8ir z 



E B ab a b (M a + M b ) + c a a t (A - M a ) 

b 



(A.4) 
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neglecting the Yukawa coupling other than top quark component a t = yf/47T 2 . Here, A 
is defined by A-term A t = Ay t . If there are no other vector-like particles below the mass 
scale of messenger field, m$, the coefficients b a , c a and B ab are given as 



and 



(&!, b 2 , & 3 ) = (33/5, 1, -3), (ci, c 2 , ca) = (26/5, 14/5, 18/5) 
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(A.5) 



(A.6) 



Then the deviation of RGE from (A.l) can be written as 

d 



dlnQ 

By using one-loop RGE for A = Aty^ 1 



]T ^a b M b + ^a t A. 



b o7T 2 87T 2 



(A.7) 



J— A = * atA + £ ^ aa M a + 0{a 2 ) 
amQ 71 a ATT 



(A. 



and one-loop part of (|A.4|), two-loop RGE (|A.7|) can be rewritten as 

'M a 



d 



— £ ^dM b + ^( 
Arc b b b 24tt V 



b 2b b 



0(a 2 ) 



where 



{d 1 ,d 2 ,d 3 ) = (26/15,6,32/3). 



(A.9) 



(A.10) 



Therefore, integrating (A.9) from messenger mass scale m$ to mz, we can obtain the gaug- 
ino mass formula at the scale mz including the two-loop effect. Because A(mz) — A t (mz), 
A(m$) — and M a (mz) — M a a~ l {m<$>) ■ a a (mz) at one-loop order, this is expressed as 



M a a-\m z ) = M a a,- l {m*) + -^-A t (m z ) 

Z47T 



4vr V V h 



Cgd b 

12h 



M b a b \m$)Aa b , 



(A.ll) 



where Aa a = a a (m$) — a a (mz)- If we write A$ = F<j,/m$ and r t = — 47rA t (mz)/A$, we 
finally obtain a following result. 



M a a-\m z ) 



C a ^ n b ( B ab 



Cad b \ > 

12X )Aat 



A$ 



(A.12) 



Since the order of r t and Aa a are 0(0.1) at most, two-loop contribution cannot become 
so large in the typical case, although it can becomes important if n a <C 1. 
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